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■ Spin relaxation in a quantum Hall ferromagnet, where filling is = 1,1/3, 1/5,..., 

, can be considered in terms of spin wave annihilation/creation processes. Hyperfine 

coupling with the nuclei of the GaAs matrix provides spin non-conservation in the 
two-dimensional electron gas and determines spin relaxation in the quantum Hall 
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c/3 . system. This mechanism competes with spin-orbit coupling channels of spin-wave 
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decay and can even dominate in a low-temperature regime where T is much smaller 
than the Zeeman gap. In this case the spin-wave relaxation process occurs non- 
exponentially with time and does not depend on the temperature. The competition of 
different relaxation channels results in crossovers in the dominant mechanism, leading 
to non-monotonic behavior of the characteristic relaxation time with the magnetic 
field. We predict that the relaxation times should reach maxima at i? ~18T in the 
i/ = l Quantum Hall system and at S ~12T for that of z^ = l/3. We estimate these 
times as ~10 — 30 /is and ~2 — 5 /us, respectively. 
PACS numbers 73.43.Lp, 78.67.De, 73.21.Fg 



2 



I. INTRODUCTION 



The two-dimensional (2D) electron gas has been intensively studied for several decades. 
The interest is stimulated by the clear manifestations of strong electron correlations in- 
cluding quantum phase transitions (like Wigner crystallization-) and, in the presence of a 
strong perpendicular magnetic field, to various features in electron transport gathered under 
the name "fractional quantum Hall effect" Transport phenomena, although of paramount 
significance for applications, provide only indirect information on such fundamental charac- 
teristics as quantum states and the energy spectrum, where optical techniques give much 
more immediate information. In particular, Raman scattering, starting from the pioneering 
works of A. Pinczuk et al,- has been successfully used to study collective excitations in two- 
dimensional electron gas created in semiconductor heterostructures and quantum wells (see 
also Ref. \4 and references therein). The position and intensities of corresponding Raman or 
luminescence lines can yield information on the energy and oscillator strengths of the excited 
states. Meanwhile an important characteristic of such excitations is also the life-time. This 
may be estimated by observation of the resonance line widths: for example the spin wave 
life-time was deduced from the observed width of the ESR resonance lines.-. Microwave and 
optical linewidths are not, however, directly related to the life-time, and usually provide 
only a very rough lower bound for this quantity. In consequence one is forced to use com- 
bined experimental methods including a time- resolved technique (see, e.g., Ref. |6|). Despite 
these experimental difficulties, growing interest in the problem of excitation life-times in a 
two-dimensional electron gas has been observed in recent years. One should mention, for 
example, recent experimental works on the observation of the spin relaxation in a polarized 
two-dimensional electron gas based on the Kerr rotation effect .- 

We study in this work the so-called quantum Hall ferromagnet where all two-dimensional 
electron gas electrons of the upper, partially filled Landau level, are in the ground state, 
with spins aligned along the magnetic field. This state obviously arises at odd integer 
fillings: u = 1,3,....- In addition, experiments and semi-phenomenological theories show 
that at some fractional fillings, namely at z/ = 1/3,1/5,..., electrons in the ground state 
occupy only one spin sub-level, and thereby the fractional quantum Hall ferromagnet state 
is also realized.-"— The quantum Hall ferromagnet possesses a macroscopically large spin S 
oriented in the direction of the field B due to negative (/-factor in GaAs structures. The 
spin wave in the quantum Hall ferromagnet may be defined as a purely electronic collective 
excitation within the Landau level which corresponds to a change of the spin numbers by 
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one: 

5S = 5S, = -1, (1.1) 

and does not alter the spin orientation of the system. (Another possible excitation in the 
quantum Hall ferromagnet is a Goldstone mode representing a deviation of S from the B 
direction which does not change the S number;— then the microscopic excitation would be a 
"zero spin exciton" corresponding to the spin change 6Sz = — 1, but 6S = 0.) This spin wave 
is also called the spin exciton, because this excitation promotes an electron to another spin 
sub-level of the same Landau level and thus an effective hole appears in the initial sub-level. 
Every spin exciton possesses energy-i^^ 

E. = ez + Sg, (1.2) 

where ez = \g\fiBB is the Zeeman gap {g ~ —0.44 in a GaAs structure); £q is the spin 
exciton Coulomb correlation energy depending on the 2D wave vector modulus q. For the 
rest of the paper it will be sufficient to consider only long wave excitations, g ^ 1/Ib {J-b is 
the magnetic length), for which the spectrum is quadratic: 

£q ^ q%/2M^. (1.3) 

Here the spin-exciton mass has the dimensionality of inverse energy.-ii^ This quantity 
has recently been measured experimentally, for u = l^^-J^ and = 1/3 fillings.— 

If there are an excessive number of spin excitons compared with equilibrium, then the spin 
relaxation reduces to an elementary process of spin exciton annihilation. The spin numbers 
are changed in accordance with Eq. (11. ip . where the energy of the annihilated excitation can 
be transferred to the emitted acoustic phonon or to another exciton due to the spin-exciton 
- spin-exciton scattering. Any spin exciton relaxation channel is thus determined by two 
necessary conditions: by the availability of an interaction that does not conserve the spin of 
the electron gas, and by a mechanism of energy dissipation making the relaxation process 
irreversible. Until now, spin-orbit coup 



conservation (see, e.g., Refs. 
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ing has been assumed to be the cause of the spin non- 
and the works cited therein). Indeed these spin-orbit 
relaxation channels are certainly dominant under the usual experimental conditions, where 
T ~ IK and 1 < i? < lOT. The corresponding calculations are in satisfactory agreement 
with the available experimental data. Here we shall extend the study of spin relaxation 
channels to include spin non-conservation by the hyperfine coupling to nuclei of the GaAs 
matrix. This has been considered previously only for the case of the Goldstone mode 
g = 0; here we consider non-zero, but small q. Our analysis shows that one mechanism in 
particular, relating to the spin-exciton - spin-exciton scattering process, should be taken into 
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account, if T < 0.1 K and magnetic fields -B > 10 T.^^ (Specifically, the necessary condition is 
T ^ ez.) To see this clearly we will analyse the spin-orbit relaxation channels - two of them 
can compete with the hyperfine coupling relaxation in the same region of temperature and 
magnetic fields. 

It should also be noted that the spin relaxation processes proceed much more slowly than 
other two-dimensional electron gas plasma relaxations unrelated to a spin change. This 
means that in any case an elementary spin exciton annihilation/creation process may be 
studied as a transition (induced by a perturbation) from an initial eigen state \i) to a 
final eigenstate |/); i.e. the hyperfine coupling relaxation mechanisms are governed, like the 
spin-orbit coupling relaxation,— li^i^"— by the Fermi golden rule probability 



where Aiji is a relevant matrix element . 

In principle, the hyperfine coupling effects are weak. The spin-orbit coupling and the 
hyperfine coupling both have relativistic origins: the former is of the first order, but the 
latter represents the second order relativistic correction to the Hamiltonian. However, the 
hyperfine interaction has some essential properties different from those of the spin-orbit 
coupling. These substantially change kinematic conditions of the spin exciton scattering 
and the dissipation mechanisms where one of spin excitons annihilates. We shall see that (i) 
first, the hyperfine coupling does not conserve total momentum of the electron system, and 
this feature leads to extension of the phase volume for the spin-exciton - spin-exciton and 
spin-exciton - phonon scatterings; (ii) second, the spin-flip process governed be the hyperfine 
interaction does not require a virtual promotion of an electron to another Landau level (this 
promotion with simultaneous spin-flip is a characteristic feature of the spin-orbit coupling 
and means a virtual conversion of the spin exciton into the cyclotron magnetoplasmon) . As 
a result, a new annihilation channel of the spin exciton scattering appears: two spin excitons 
can be scattered by each other, within the same Landau level, directly due to the hyperfine 
interaction - finally one gets a single-spin exciton state possessing the combined energy. 
This kind of scattering, as in the case of scattering caused by disorder,— is kinematic: 
the transition matrix element does not contain the Coulomb constant between bra- and ket- 
vectors - the scattering is possible because the spin excitons are not actually elementary 
Bose particles but possess an internal degree of freedom and thus have a "memory" of the 
Pauli principle for the primary electron system. Thus in spite of small hyperfine coupling 
constant, the hyperfine coupling channel competes with the spin-orbit ones and can even 
dominate. 



Wfi = {2n/h)\Mfi\H{Ej-Ei), 
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The next section of the paper is devoted to formal description of the system where we 
present the Hamiltonian and the basis of exciton states (excitonic representation). In section 
III we study the hyperfine coupling relaxation mechanisms when the spin exciton annihi- 
lation/creation is determined by the spin-exciton - spin-exciton scattering including the 
dynamic and kinematic scattering channels. For this process the relaxation rate is propor- 
tional to the spin-exciton number squared, and therefore the relaxation is non-exponential 
with time. (In principle, it becomes exponential when the spin exciton number approaches 
its equilibrium value, but the final exponential stage cannot, in fact, be observed under the 
condition T ^ ez.) We discuss also in Sec. Ill possible relaxation processes, related to the 
hyperfine coupling and phonon emission/absorption, comparing them with other relaxation 
mechanisms. Section IV is devoted to the spin-orbit relaxation channels relevant to the con- 
sidered region of strong magnetic fields and low temperatures. These spin-orbit mechanisms 
are also related to the spin-exciton - spin-exciton scattering but determined by two different 
dissipation processes - via coupling to a smooth random potential or by coupling to phonons. 



In Sec. V we discuss the results of our study. The main result consists of the interplay 
of different relaxation processes. We compare those due to the hyperfine coupling and spin- 
orbit interactions and, summing all relaxation channels, calculate the total characteristic 
inverse time. In this "interplay regime", where the spin-exciton - spin-exciton channels 
dominate, the relaxation occurs non-exponentially, and the effective relaxation time reaches 
its maximum ~ 1 — 5 /is depending on the Landau level fillings. Nevertheless, the relevant 
region of parameters T and B is not too extreme and experimentally quite accessible. 



Note that in this paper we do not study the situation where the Goldstone condensate 
of "zero spin-excitons" arises,—!^, i. e. where there would be a rotation of the direction but 
not a reduction in the magnitude of the total spin of the system. Here we consider instead 
relaxation where there are, at low temperatures a bulk number of spin excitons arising from 
an intensive external (e.g., optical) excitation. The initial state at low temperatures should 
be described as a (metastable) "thermodynamic condensate" of spin waves with non-zero, 
but small wave-vectors limited by the uncertainty determined by disorder. We think that 
such situation, where most of the spin exciton annihilation/creation events happen within 
the thermodynamic condensate, is realisable experimentally. 
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II. FORMAL STATEMENT OF THE PROBLEM: THE HAMILTONIAN AND 

THE BASIS OF EXCITON STATES 

Our system consists of two components: electrons belonging to the two-dimensional elec- 
tron gas and nuclei of Ga and As atoms. In addition, we consider piezo- and deformation 
couplings of the 2D electron gas electrons to the lattice, which are reduced to electron-phonon 
interaction. So, the Hamiltonian used is as follows: 

H = H, + ezSz + Ant + E + E ^i-Ph • (2.1) 

j 3 

Here = + H^o is single electron Hamiltonian including the spin-orbit cou- 

pling part (q = —iV + eA/cK); Sz = ^^a^J^ /2 ; subscript j labels electrons. The third 
term describes Coulomb energy of the e-e interaction, the fourth is the hyperfine interaction 
of electrons with nuclei, and the fifth is the operator of electron-phonon interaction. If one 
holds Hso = 0, we can omit the orbital single electron energy terms — all states relevant to 
our problem belong to the same Landau level, and therefore have the same orbital energy 
equal to huJcuN^ {uc is the cyclotron frequency, A/^ is the Landau level degeneracy). We 
ignore also the energy of nuclei which consists of the contribution due to their interaction 
independent of the electrons, and of the nuclear Zeeman energy. Variations of both, associ- 
ated with change of nuclear spins, are negligibly small owing to the tiny nuclear magnetic 
momentum. 

In the following three subsections (A, B, and C) we neglect the spin-orbit coupling. The 
spin-orbit Hamiltonian and spin-orbit corrections, written in terms of the representation 
used, will be given in the subsection D. 



A. Electron system. Excitonic representation 

We now present the basis set of states diagonalizing the first three terms of the Hamilto- 
nian |2II]to leading order in parameter = [ae^ j kIb) / h^c considered to be small (a < 1 is the 
averaged form-factor arising due to finiteness of the 2D layer thickness; k is the dielectric 
constant). We do this by analogy with previous works,— "i^i^S^— defining the spin exciton 
creation operator— 

2-. = ;;;^E--"^t.f%-f ' (2-2) 

where Op and hp are the Fermi annihilation operators corresponding to electron states on the 
upper Landau level with spin up (0=^) and spin down {h=\), respectively. Index p marks 
intrinsic Landau level states which have wave functions ?/'„p(r) = (27rA/'^)~"'^/'^e*^^V9„(p+x) in 
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the Landau gauge. [^Pn{x) is the oscillator function, where n is number of the upper partially 
filled Landau level; in Eq. fl2.2p and everywhere below we measure length in the Ib units 
wave vectors in the 1/Ib ones.] In the odd-integer quantum Hall regime, operator (12. 2p acting 
on the ground state yields the eigenstate of the first two terms of Eq. (12. ip . namely: 

[6z5. + ^int, QLj|0) = (6z+i:,)QLq|0), (2.3) 



where |0) = 1 1) t? ••• t)- This basic property of the exciton state, Q)jf,q|0), is asymptotically 



exact to first order in r^.. After the introduction of intra-sublevel operators A\^ = N'^ 2Lq 
and B]^ = A/'^ ^^^SLq' obtain a closed Lie algebra for these exciton operators.—"^— The 
commutation identities needed in our case are 

\0 0+ 1 — P*('li^'l2)z/2 /t _ p-j(qiXq2)z/2K? 

[isiqi, isiqaj — C ■^qi-q2 '-'qi-q2' (0 A\ 

,~.(q.Xq.)./2[^t^^ Qt J = _e^(..Xq.)./2^^t^^ Qt J = -U.'QI^^ 

(Here and below we omit the subscript ah at the Q-operators.) Note that the commutation 
algebra (12. 4p is neither purely Fermionic nor Bosonic. 

The interaction Hamiltonian H\y^t = |/ 'iri(ir2 ^'^(r2)\l/^(ri)f/(ri — r2)^(ri)\l/(r2) may be 
expressed in terms of the exciton operators.—*^ If we keep in i^int only the terms relevant 
to our problem, it takes a very simple form 

^int = ^ E ^(^) + 2^iSq + ^i^q) , (2-5) 

q 

Here W{q) = U{q)[f{q)]^, where / = e"'?'/^ if z/ < 1, or / = e-'''/^[Lfe(gV2)] if u = 2n+l (L„ 
is the Laguerre polynomial). U{q) is the Fourier component of the 2D Coulomb interaction 
function: U{q) = {e^ / nlBq) ff (izi(iz2e~'^'^^^^^'|x(-2i)Plx(2^2)P, where x(^) is the dimensionally 
quantized wave-function of an electron sized in the 2;-direction. 

In contrast to integer quantum Hall ferromagnet, the use of the excitonic basis Qq|0) 
presents only a model approach in the case of fractional quantum Hall regime. Generally, at 
fractional filling, spin-flip excitations within the same Landau level might have many-particle 
rather than two-particle nature because the same change in the spin numbers (II. ip may be 
achieved with participation of arbitrary number of intra-spin-sublevel excitations (charge- 
density waves). These waves are generated by the operator Al^ acting on the ground state 

|0) = I 'I', .. t, .. 'I' ).- The result is trivial in the case of integer u (^q|0) = 5q,o|0)); however, 
states of the Ql^^Al^^Al^^...\0) type might constitute a basis set if one studies the spin-flip 
process at fractional u. On the other hand, a comprehensive phenomenological analysis-*^ 
suggests that even the spin-flip basis reduced to single-mode (single-exciton) states would 
be quite appropriate, at least for the lowest-energy excitations in the case of fractional 
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quantum Hall ferromagnet. This single-mode approach is indirectly substantiated by the 
fact that the charge-density wave has a Coulomb gap ^ which is substantially larger than the 
Zeeman gap ez- Hence for a fractional quantum Hall ferromagnet, just as in Ref. 0, we will 
consider the simple state Qq|0) to describe the spin- flip excitation. However, the calculation 
of (0|^q^q,|0) is required for the following. Now this expectation is not simply equal to 
5q o5q',0) but is cxpresscd in terms of the two-particle correlation function g{r) calculated for 
the ground state: 

{0\A^A[,\0) = [27rz/^(g)e«'/2^l] V,q < 1) • (2-6) 

Here 'g{q) = j^^^ j g{r)e~'^'^^(P'r is the Fourier component. The function g{r) is known, 
e.g., in the case of Laughlin's state.-ii^ If the ground state is described by the Hartree- 
Fock approximation, we have simply 2Tcg= [N',j,5q^Q — e'~'^^^'^^, which does not depend on v. 
Furthermore, at odd-integer filling factors this Hartree-Fock expression becomes a Fourier 
component of the exact correlation function. In the latter case one should also formally set 
z/ = 1 in Eq. (12. 6p Note also that the exact equation (0|^Jj|0) = z/'(5qo holds, where we set 
v' = v ii V <1 but z/' = 1 if z/ = 3, 5, ... 

With the help of Eqs. (12. 4p and (12. 5p one can check Eq. (12. 3p in the case of odd-integer 
V. If V is fractional, the Coulomb exciton energy within the single mode approximation 
is defined as £q = {0\Qci[Hint, Qq] |0)/(0| QqQ]j|0)."'^° As a result in both cases of integer or 
fractional u < 1 one obtains for small q the quadratic dispersion law (II. 3p with the spin 
exciton mass 

= IJ^^ W{q)q' (^1 - ^(Ol^q^t |o)^ dq . (2.7) 

We have employed the rule for change from summation to integration over the 2D vector q: 
Y.^...=M^j ...qdqdct>/2'K. 

B. Nuclear system. Hyperfine coupling in the excitonic representation. 
Hybridization of spin-exciton and nuclear spin-flip states. 

The general expression of the hyperfine coupling Hamiltonian is simplified in the case 
of interaction with nuclei in a semiconductor matrix. ^^"^^ As this simplification is valid in 
the 2D channel of a quantum well we may directly start from the well known expression for 
contact interactions of electrons with nuclei 



i^hf = ^5^A„vI/*(i2„) (/^"^(t) ^>{R^) , (2.8) 



(see for example Ref. |30| and references therein.) where I and Rn are spin and position of 



the n-th nucleus and '^{R) is the envelope function of electron [R = (r, z) is the 3D vector]. 
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Both Ga and As nuclei have the same total spin: I^^ = I^^ = 3/2. In Eq. fl2.8p Vq is volume 
of the unit cell. The parameter An, being inversely proportional to Vq, really depends only 
on position of the Ga/As nucleus within the unit cell. For the final calculation we need the 
sum Aq^+A\^. If t>o is volume of the two atom unit cell then Aq^+A\^c::^4:- 10~^ meV^. (See 
Appendix A.) 

- (n) - - - 

Now we rewrite / -ct as /^a^ + /+(3"_ + /_(3'+. Then omitting the J^a^ term due to its 
irrelevance to any spin-flip process and substituting in Eq. f l2.8p the Schrodinger operators 
^^{R)=xiz) Y,p (o-l + bl) V'p(r) and ^(i?)= (^^)^ instead of and ^, we come to 

Hm = ^Y1 ^J^^Pi E \xiZn)\'rr,Mn, Yn)i'pAXn, Yn)Aj^^^ + H.c. (2.9) 

p\,P2 n 



Substitution of the equation 



igx{p2-<ly/2) 



bW = E 7t^K,p^p^ Qi ' (2-10) 

q V^^ 

which is simply inverse to Eq. (12. 2p . yields after summation over pi and p2 the hyperfine 
coupling Hamiltonian in the excitonic representation: 

Hm = ^^/°^ E/(g)Q-^E-^"l^(^-)l'^^'''"-^-^ + H.c. (2.11) 

A set of the spin numbers {M} = (Mi, M2, ...M^j, ...), where every M„ may take on 
values —3/2,-1/2,1/2,3/2, completely determines the state of the nuclear system. The 
state where 2D electrons are in the ground state and nuclei in the state {M} we denote as 
|{M},0). By applying the lowering/raising operator I^'^ to this state, we obtain 



4"^|{M},0) = ^QtM„) Q±M„)|{M};^,0>, (2.12) 

where {M}J = (Mi,M2, ...M„=pl, ...). Let us find the hyperfine coupling correction to the 
normalized spin exciton state Q^\{M},0) / y/iT'. Considering operator (12. lip as a perturba- 
tion we obtain with the help of Eqs. (12. 4p and (I2.12p : 



|SE,q) = Qi|{M},0)/vV 

(2.13) 



4 
in 



^|||^E-^l.(-.0|V-\/(|-M„)(|.M,,)|,M,-,0> 



In the same way we find the corrected nuclear 'spin-turned' (NST) state 
|NST,n) = |{A'/};,0) 



(2.14) 
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where we consider M„ > —3/2, and use notation {M}~^, = (Mi, ...M„ — 1, ...M„/ + 1, ...) 
meaning by that {M}~+ = {M}. [The n' = n term in the sum of Eq. fl2.14p contributes to 
the transition matrix element relevant to some spin exciton relaxation processes.] 

The hybridized states f l2.13p and f l2.14p diagonalize the first three terms of the Hamiltonian 
(12. ip to the first order in hyperfine coupling. Correspondingly, these have energies E^{q) and 
(counted from the energy of the |{M}, 0) state) within the approximation neglecting energy 
corrections of the second order in hyperfine coupling and small magnetic energy corrections 
related to changes of M„ momenta. 



C. Electron-phonon interaction in the exciton representation. 

The Hamiltonian of the interaction of electrons with 3D acoustic phonons is written as:— 

^e-ph = ^'^^^^ ^k,s'He-ph(q) + H. c. , (2.15) 

where L"^ = 2tiM^I\ is the 2D area, and is the dimension of the sample along z, 

•He-ph(q) = j e^'''-^t(r)^(r) d\ , k = (q, ; (2.16) 

Pk,s is the phonon annihilation operator (index s denotes possible phonon polarizations: 
the longitudinal / or one of the two transverse polarizations ti or ^2); and ?7's(k) is the 
renormalized vertex which includes the fields of deformation (DA) and piezoelectric (PA) 
interactions. The integration with respect to z has been already performed and leads to the 
renormalization t/'s(k) = f/s(k) J x*{z)^^^''^x{z) dz 

The isotropic model for the phonon field— enables us to take into account the deformation 
and piezoelectric couplings independently. We further use the approximation where we take 
no difference between longitudinal and transverse sound velocities. For the three-dimensional 
(3D) vertex one needs only the expressions for the squares,— 1^ 

|f// = 7r£ph(k)/p3r,(k), (2.17) 



where the phonon energy is = hc^/kf+(f /Ib (we recall that K, and q are dimensionless) , 
Pq = 2.52 ■ lO^cm"^ is the material parameter of GaAs (see Ref. l33l ). The longitudinal Ti(k) 
and transverse rf(k) times are the 3D acoustic phonon life-times (see Appendix A). These 
quantities are expressed in terms of nominal times and Tp characterizing respectively DA 
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and PA phonon scattering in three-dimensional GaAs crystal. (See Appendix A, Ref. 
and cf. Ref. i33 .) 

The dimensionless operator "He-ph in terms of the excitonic representation has the follow- 
ing simple form (cf. Ref. 



11 



n^-pM = f{q)M^ (A + ^q) • (2-18) 

D. The spin-orbit coupling in the excitonic representation. 

If considering the spin-orbit coupUng, we will ignore the hyperfine coupling but take into 
account the Hso operator in the single electron part Hi of the Hamiltonian (12.11) : 

-f^so = a (q X (t)^+/3 (gj^(Ty-ga;(Ta;) , q=-iV + eA/c^. (2.19) 

This operator, specified for the (001) GaAs plane, represents a combination of the Rashba 
term (~ a) and the crystalline anisotropy term (~ P)— and does not violate translational 
symmetry.— 

Now it is convenient to use a bare single-electron basis diagonalizing Hamiltonian 
q^/2m* + Hso- To within the leading order in the ifso terms we obtain 

= ( ^ , *" , ^ , ) and = (""V^^'-P + "V^-A-..) . (2.20) 

where u and v are small dimensionless parameters: u = (3\/2/IbItiOc and v = a 
Thus the single-electron states acquire a chirality a or 6 instead of spin quantum number, 
and the spin flip corresponds to the a — > 6 process. The definition of the spin exciton creation 
operator formally remains the same [Eq. (12. 2p ]. however the and bp operators describe 
annihilation in the states (I2.20p now. 

When being presented in terms of basis states (I2.20p . spin operators J \E'"'"S^\E'(i^r and 
/ '^'^Sz'^cPr [where = J2p{(^p'^pa + bp'^pb)] are invariant up to the second order of u and v. 
However, in the excitonic representation the interaction Hamiltonian i^int and the electron- 
phonon coupling operator acquire terms proportional to u and v, which are additional to Eqs. 
(12. 5p and (I2.18P respectively.—"—!^ These terms correspond to creation and annihilation of 
spin excitons in the system: 

HL = - ^^~W{q) K + K) 2q + H- c, (2.21) 

q 

and 

?^Lph(q) = K"f{q) - vg-) + h. c. (2.22) 

[q^ = Tiiqx±iqy)/V^- 

We can also take into account the presence of an external smooth random potential ^p{r). 
This is assumed to be Gaussian and defined by a correlator K{r) = {(p{r)ip{0)). By choosing 
{^p{r)) = 0, , the correlator is K{r) = exp (— r^/A^), in terms of the correlation length 
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A and the amphtude A. The smooth random potential can act as the rate-hmiting process 
in the energy dissipation which makes the spin-flip process irreversible. (p{t) formally is 
analogous to frozen field of phonons having zero frequency. A static potential cannot cause 
dissipation alone: physically the random potential (mixed with the spin-orbit term) causes 
spin-flip and breaks momentum conservation. The actual dissipation comes from other 
interactions that do not change the spin: electron-electron and electron-phonon interactions 
that occur on a faster time scale and render the process irreversible. Therefore, using again 
the Eq. (12.201) basis set and Eq. (12.101) . we obtain the operator in terms of the excitonic 
representation. The part responsible for a spin-flip is— 

^' = E (^«?+ - ^^-) 2q + H-c. , (2.23) 

q 

where ^ is the Fourier component [f{r) = X]q^(^)^"''']- 

III. THE SPIN-EXCITON - SPIN-EXCITON SCATTERING RELAXATION 
CHANNELS GOVERNED BY THE HYPERFINE COUPLING 

The 6Sz = —l hybridized states (I2.13P and (I2.14p diagonalize the Hamiltonian ^int+-^hf; 
but the 6Sz = —2 states Q^JSE, and Q^INST,?!) do not. (Here by we mean the 
total spin number of the combined nuclear and electron system.) The problem may be 
formulated in terms of a scattering where the double exciton state Qj^JSE, transforms 
to the single exciton one Q^|NST,n). Since the hyperfine coupling energy is neglected, the 
energy conservation law takes the form 

E,{q,) + E^iq^) = EM ■ (3-1) 

It determines the modulus of the spin exciton momentum q in the final state, and in partic- 
ular means that q cannot be equal to qi or q2. 

A. Kinematic scattering 

The transition matrix element Mif in Eq. (II. 4p has to be found to first order in the 
hyperfine coupling. Therefore in the case of the kinematic scattering, where Aiif represents 
an expectation value (bra|ifhf |ket) calculated directly for the hyperfine coupling operator, 
the ket- and bra- vectors are determined only by the main components of the Jgg^ and 
Q^INST,?!) states without any hyperfine coupling corrections. Namely, taking into account 
that the initial double-exciton state and the final single-exciton one have to be normalized. 
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we should calculate the kinematic scattering matrix element 

-MS"(qi,q2,q,ri) = {Q,{M}-\Q^H^,Q\^Q\^\{M},Q) /v"'\ (3.2) 

After substitution of Eqs. (12.111) and (12.121) . this is reduced to calculation of the four- 
operator expectation value (C.l) (see Appendix C). Note that were the Q-operators usual 
Bose operators, the expectation (C.l) would simply be equal to 5q',qi5q,q2+ 5q',q2^q,qi) the 
conservation condition (13. ip could not be satisfied and therefore the kinematic scattering 
channel would not exist. Therefore only due to the non-Bose nature of the spin-exciton 
states does this relaxation mechanism take place. 

We should keep in A^^™ only the main terms contributing to the final result, namely to 
the relaxation rate calculated on the basis of the Fermi golden rule (ll.4p and subsequent 
summation over the qi, q2 and q statistical distributions. These are terms to the lowest 
power of qi q2 and q. They give the exact result to leading order in the small parameter TM^. 
(T is the temperature, characteristic values of the momenta are gi, g2, q ~ y/TM^ ^ 1.) In 
particular, one finds that the ~ v/M^ terms in Eq. (C.l) give the strongest contribution, 
and the ~ (0|^.,,^^ |0) terms may be neglected.— In addition, the terms where q = qi or 
q = q2, are omitted due to the 'selection rule' determined by Eq. (13. ip . As a result we obtain 



A,-(,„,„,.„) /(!±^pzM,.<.«..«, . (3,3) 



B 



B. Dynamic scattering 

If studying the dynamic scattering, one should take into account that the Coulomb in- 
teraction operator (12. Sp . acting on a certain state, does not change the number of the spin 
exciton operators determining this state, - i.e. this number must be the same in the bra- and 
ket-states contributing to M.if = (bra|i^int|ket). Furthermore the Coulomb interaction does 
not change the total momentum of the electron gas, - it too must be the same in the bra- 
and ket-states. Therefore, again only the ~ Qq^ Qqa |{^}) 0) component should be kept in 
the initial state Qj^JSE, q2). [The hyperfine coupling correction component can contribute 
only to the transition where q = qi, which is forbidden due to Eq. (13. ip . 

The single exciton state Qq|0) diagonalizes the Hamiltonian -ffint, but the double exciton 
state QjjjQjjjlO) does not. The latter is in fact an 'almost' eigenstate. Indeed, even at 
odd-integer v we have [cf. Eq. (12. 3p 

[i^int, QlQlm = {£,. + £,^ QlQ\M + [[^int, Qll QiJ |o), (3.4) 

where the double- commutation term arises due to the interaction between the spin excitons. 
It can be routinely calculated with the help of Eqs. (12. 4p and (12. 5p . see Eq. (C.2) in 
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Appendix C. The norm of this term and the averaged spin-exciton - spin-exciton interaction 
energy (0|Qq2QqJ [Hint, Q\^],Q\^ |0), both vanishing if gi = or = 0, are respectively 
< {ae^ / hIb)-M^ ^^"^ and < [ae^ j kIb)I if The latter estimation quite corresponds 

to an effective mean dipole-dipole interaction of two spin excitons sized within the area 
2t[1\M^. [We recall that each magneto-exciton possesses a dipole momentum equal to e(q x 
z)l\ (in common units).—] 

It follows from the above that for the dynamic scattering process we choose the ket- and 
bra-vectors as 



|ket) = Q^^Q^J{M},0)/z/', 
and (bra| = (n, q| 



^^0 



(3.5) 



implying that only the hyperfine coupling correction term is relevant in the final normalized 
state Qj^|NST, n)/z/'"^^^. The matrix element meant to be calculated is 



Kr(qi,q2,q,^) = (^,q [[Hint, Ql^l Ql\ {m},o)/u'. (3.6) 

By analogy with the kinematic scattering, we keep in A^fj° only terms to the lowest power 
of qi q2 and q. Using sequentially Eqs. (C.2), (C.l), (12. 7p and (3.1), we find 

Mi} (qi, q2, q, n) = — -— . (3.7) 

+ qiq2 - q(qi+ q2) 



C. The relaxation rate 

To calculate the spin-wave relaxation rate, one should know the distribution A^q of spin 
excitons over the q wave numbers. Although exciton operators (12. 2p are non-bosonic, the 
spin exciton obey Bose statistics, because their number in any state determined by a certain q 
may, in principle, be macroscopically large. At any moment the spin excitons distribution is 
in quasi-equilibrium and characterized by a chemical potential fi<ez- (The thermodynamic 
equilibrium is established much faster than spin-flip processes occur.) Initially the total 
number of spin excitons N':^ = v'M^/2—S is actually determined by a short external optical 
impulse, and its value might be even more than the critical value 

A4c=Ar, r .f%^_ =U,M^T \ql/2M^T - In (e^oV2M.T_^^j ^ ^3 

J qo 6Xp(^C.q/ij i L V /J 

where we have used the quadratic approximation (11. 3p and designated as go a lowest limit 
of possible nonzero values of g. Any violation of the translation symmetry contributes to 
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the estimation of go- For example in the ideaUy clean case qo ~ ^/L, where L ~ ^/J^ is 
the linear dimension of the 2D system. A more realistic estimation can be made if one takes 
into account the presence of a smooth random potential, then go ~ Afx^sA/A, where A is 
the potential amplitude (A^l/Mx), and A is the correlation length (A^Z^)-— In practice 
go < 0.01. If A/'x > A/'xc, then the bulk number of spin excitons with nonzero but momenta 
|q|<go form a thermodynamic condensate. The specific g-distribution of excitons within 
the condensate plays no role; however we may write 

r{0} 



if q G {0} 

l/[exp{SjT)-l], if q^{0} 



(3.9) 



(q G {0} means belonging to the thermodynamic condensate). The number of the conden- 
sate excitons is thus A4— A/'xc = X] ^q""*^- During the spin exciton relaxation process the 

qe{0} 

condensate is depleted, and when A/'x < A/'xc we have: 



Arq = l/[exp {£^ + ez-fx)/T-l] 



(3.10) 



with chemical potential equal to 

/(i = ez + Tin 



exp 



A4 



(3.11) 



AT^MxT^ 

[In the vicinity of ez the value /i is determined with an accuracy: — ez| ^ min(gQ/2Mx, T).] 

The = equation determines the equilibrium spin exciton number: ATio) = -AT^MxTln (1 - e-^^/^). 

The spin wave relaxation rate is defined as the difference between the fluxes of annihilating 
and created spin excitons in the phase space: 

"^■^^ ^ J2 Si^u^2)[N^,N^,a + Nu) - N^2il + N^,){1 + N^,)], (3.12) 



dt 



qi,q2 



where Ni2 = l/ [e^^n+^^2+^z)/T if A4>A4c or A^ia = 1/ [e^^'i+^^a+^^^-M)/^ - 1] if A4<A4, 

and the summation over final state values q is performed by calculating 

27r I ^ 

5(qi,q2) = -r-^^ -M^7(qi,q2,q,r^) + A^f/°(qi,q2,q,ri) 5{£g- £q,- £q^- ez) 



n q 

t;2Mx [1 + J^iql ql 0, 2Mxez)] / 3 



where cj) is the angle between qi and q2. 



-MjA^|x(^n)r (3-13) 



2 '7 V2 
[1 + J^{ql ql 0, 2Mxex)] /AT^Tm , 



xy cos^0 {/3 + X + y) 



+ f3{x + y) + xy cos 
1 _ 5t;oMx(A|, + Ai 



2^1 3/2 ' 



and 



(3.14) 
(3.15) 
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[we have kept in T only terms nonvanishing after averaging over the qi and q2 directions 
when in Eq. fl3.12p ]. The summation over n in Eq. f l3.12p has been performed for the case 
of unpolarized nuclei. In addition the correlation length of the spatial nuclear momenta 
distribution has been considered to be smaller than the magnetic length Ib and conventional 
width of the two-dimensional electron gas: d=(^j \x{z)\'^dz^ ^. (This value is certainly not 
equal to the quantum well width dqy/, but constitutes a fraction of it, e.g.: d/dqy^ — 1/3.) 

The rate —dM^/dt is completely determined by Eqs. fl3.9l) - fl3.15p . In the following calcu- 
lations we use the following: (i) the kinematic and dynamic scattering fluxes simply add, as 
independent contributions to the total rate; (ii) in the case of T^ez the contribution to the 
rate due to the dynamic scattering relaxation flux is negligibly small; the same result is found 
if one of spin excitons in the initial state belongs to the thermodynamic condensate (i.e. qi 
or/and q2 G {0}); (iii) S'(qi, q2) does not depend on qi and q2 for kinematic scattering, and 
the summation in Eq. (13.120 reduces to 

[•••] = - E ^i2(l+A^q.+A^qJ • (3.16) 
qi,q2 qi,q2 

In the T>ez region the spin-orbit relaxation channels are much more intense than the 
considered hyperflne coupling channel (see the next sections), and both spin-orbit and hy- 
perflne coupling relaxation mechanisms compete with each other only in the T ^ ez case. 
Therefore we specifically study this situation. Then the dynamic spin-exciton - spin-exciton 
scattering is neglected, and the spin exciton creation term in Eq. (I3.16P may be presented 
as - Eq„q. ^i2(-) ^ ~e-^^^ iV^,iV^.(l+iVq,+iV,,J/(l+iVqJ(l+iVqJ. In the » T 

case this term is a negligible quantity compared to A/^. If we consider /i<T, then the term 
is equal to — A/'xA/'i*^'' . So, if T <^ ez, then for any relation between /i and T one finds that 
Eq. (13.120 reduces to 

- dnjdt = n^[n^ - 4°)] /2tm (T<ez) (3.17) 

[nx(t) = M^{t)/M4, and ni°'' = N'^p^/Af^ to note the spin exciton concentrations]. In fact 
under the conditions considered, the observable relaxation process is completed while still 
ni°^t/2rhf<l, then 

n,(t) = ^^i^l- . (3.18) 

l + nx(0)t/2rM ^ ' 

This law is independent of the temperature but depends on the magnitude of the initial spin 
excitation nx{0). The effective relaxation rate is ~ nx(0)/2rhf < 0.1/rhf (if one assumes that 
nx(0)<0.1). 
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D. Spin exciton relaxation due to hyperfine coupling together with the interaction 

of spin excitons with acoustic phonons 

In principle, the spin-exciton - phonon coupling mechanism participates both in the spin- 
exciton - spin-exciton annihilation scattering and in the single-spin exciton one. However in 
the case of spin-exciton - spin-exciton scattering this relaxation channel represents only a 
small correction to those studied in the previous subsections, proportional to electron-lattice 
coupling constants. Let us estimate the spin exciton-phonon relaxation governed by the 
single-exciton annihilation mechanism. We need to calculate the transition matrix element 
■Mx-ph between the state |ket) = |SE, qi) and some of final states |bra) = Pj^ ^|NST, ra) 
for the exciton-phonon operator determined by Eqs. (12.151) and fl2.18p . Now the energy 
conservation law reads -E'x(qi) = hck/lB, where k= {kz,q). Meanwhile the q = phonons 
do not contribute to the relaxation process, because action of the 7{e-ph(0) operator fl2.18p 
on the |SE, 0) state is reduced to multiplication by a constant - hence A^x-ph = due 
orthogonality of the |SE, qi) and |NST,n) states. If g 7^0 then the contribution to Aix-ph 
is determined only by the first component of the ket-state |SE, qi), namely by commutators 
['He_ph(q), QqJ \{M},0) /y/i^. The latter according to Eq. (12.161) and commutation rules 
(12. 4p vanish in case qi = being proportional to qxqi at small qi. This issue is a key point: the 
matrix element squared |A^x-ph(qi)P is proportional not only to the small constants of the 
hyperfine coupling and electron-phonon coupling but also to the temperature (more exactly 
to the small dimensionless parameter M^T). As a result, making computations similar to 
those made above, we finally obtain a relaxation rate linear in n^: — c/Arix/rft = Arix/rhf-ph 
[A Tlx to note the difference nx— with the characteristic inverse time 

1 u'voM^jAl^ + AjJezT ^ 
Thf-ph hd%dplTD 

(under the considered conditions predominantly the deformation part of the e-phonon cou- 
pling contributes to the result). This value is much smaller than the inverse time given by 
formula (I3.15p . Much more important is comparison with another value governing also the 
single-spin exciton relaxation process related to phonon emission: namely, a certain char- 
acteristic inverse time l/xso-ph can be calculated in the case where spin non-conservation 
instead of the hyperfine coupling is determined by the spin-orbit coupling.— 1^ It is found 
that at any parameters 1 / Xso-ph is much larger than 1 / Thf-ph (by two or three orders of mag- 
nitude). We conclude that spin exciton relaxation channels appearing due to the hyperfine 
coupling together with electron-phonon coupling are very slow and may always be neglected. 
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IV. THE SPIN-ORBIT RELAXATION CHANNELS 



The spin-orbit relaxation channels, similarly to the hyperfine coupling mechanisms, may 
be subdivided into the two spin-exciton scattering channels and the single spin exciton ones. 
Among them there is a strong spin-exciton - spin-exciton scattering process actually respon- 
sible for the spin exciton relaxation under the conditions of published experimental studies,-*^ 
namely at T ~ 1 K and i? < 10 T. This is the spin-exciton - spin-exciton dynamic scattering 
where the spin- flip is determined by the transition matrix element (fin|if'jj^.|ini) calculated for 
operator fl2.2ip . and states |ini) = Q^^Q^JO)///' and |fin) = Qq|0)/-\/z^. Being constrained 
by energy E{qi)+E{q2) = E{q) and momentum conservation qi+q2 = q, this process occurs 
if qiq2 = ezMx; i.e. the phase volume of the scattered spin excitons is essentially restricted. 
In particular, if the scattering spin excitons belong to the thermodynamic condensate, this 
relaxation mechanism is switched off. In fact the dynamic relaxation channel works well only 
when T>ez, giving the relaxation time ~10ns.— "^^"^ However if T ^ ez, the characteristic 
time is drastically extended, as it is proportional to the double exponent ~ ^'^^1'^ (^see Ref. 
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). Therefore, studying exactly the T <^ ez case where the spin-orbit and hyperfine coupling 
relaxations are competing, we consider the spin-exciton - spin-exciton kinematic processes 
provide more intense relaxation. In the excitonic representation these are determined by 
operators fl2.22p and (12.231) which do not conserve the number of spin excitons. 



A. Relaxation via a smooth random potential 

The spin-orbit relaxation channel in presence of a smooth random potential, is again gov- 
erned by the kinetic equation (13.121) where 5'(qi,q2) = (27r//i) |-^i/''(qi) q2, q) | ^{^q — 
£q^~£q^—ez), where M.^^j' = (fin|(^'|ini) with initial and final states |ini) = Ql^^Ql^^\0) / u ' and 
|fin) = Qq|0)/v^ respectively. Taking into account that Sq^.Sq^ <^ ez, the argument of the 
^-function may be set £q — ez, and using Eq. (C.l), we obtain the squared matrix element 

11 2 2 

[A^-^'^l =2{u'^ + v'^) l^'^l^ where q* = y/2M^ez, and the scattering probability 

independent of qi and q2: S = l/A/^Tg^P. The characteristic inverse relaxation time is 

l/r,7 = WTr'iu^ + v'')MlezK{q*)/u'h. (4.1) 

Here K{q) stands for the Fourier component of the correlator. If the latter represents a 
Gaussian function (see Sec. II-D), then 'K{q*) = ttA^ exp {-M^ezA"^ /2ll).^^ We note that it 
depends exponentially on the magnetic field squared: r^e^'^^^ (the spin exciton mass is as- 
sumed to be independent of B). As mentioned earlier, this time is assumed much longer than 
the times of thermalization and therefore determines the relaxation while the irreversibility 
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occurs due to the fast thermalization. The relaxation rate can then be calculated as in Sec. 
III-C. Whether or not the thermodynamic condensate exists, it is governed by equation 

- dnjdt = [n^ - ] /2r,7 (4.2) 

differing from Eq. (13 . 1 7p only by the replacement of Thf with r^^P. Likewise one obtains Eq. 
(I3.18P with the same substitution. 



B. Electron-phonon coupling mechanism of the dissipation 

We study in this subsection the spin-exciton - spin-exciton scattering process, where there 
are two spin excitons in the initial state and a single spin exciton plus an emitted phonon in 
the final state. (For a discussion of single spin-exciton annihilation due to phonon emission, 
see comments at the end of Sec. III-D). In this case the conservation laws read: 

qi + q2 = q + qph and 

Now the kinetic equation for annihilated and created spin excitons is 

- ^ = ^ E ^(^1' q) [A^qi A^q. (1 + A^q + A^ph) ' iVqiVph (1 + N^, + iVqJ] , (4.4) 

qi,q2,q 

Due to the T ^ ez condition we can neglect values Sq-^^ and Sq^ in the £'(qi)+£'(q2)— -E(q)— eph 
argument of the (5-function when calculating the scattering probability, therefore 

'^(qi,q2,q) = Y Yl l-^x-ph(qi)q2,q,^^,qph,s)|^(5(ez-^g-^c-yA;2 + g2^). (4.5) 

The matrix element is A^x-ph = (fiii|-ffe-ph|iiii), where the electron-phonon Hamiltonian is 
presented by Eqs. (l2.15p -( |2rT7j) (with change from 7/c-ph to 'Hp_pjj, see Eq. (12.220 ). and 
bra- and ket- vectors are |fin) = -Pfc^ s^ql^)/^^^ 1^^^) ~ 2qi2q2|0)/^' respectively. 
Using Eq. (C.l) we keep again only terms ~ v'/M^ contributing to the result in the leading 
approximation. Finally, by doing in the spirit of manipulations above, one obtains the 
relaxation rate (14. 4p in the form 

- dnJdt = njnx- j/Sr^P'^ (4.6) 

similar to Eqs. (I3.17P and (14.21) . Now the temperature-independent constant characterizing 
the rate is — 
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where 

^(0= / dx{x-x'^)/^/l-2^x/{l-xy 
Jo 

[xq = 1+^ — a/^^ + 2^]. In the derivation we have set 1/ta ~ 1/td, because estimation shows 
that contribution of the deformation couphng is dominating the polarization one under the 
considered conditions (cf. Sec. III-D). Unhke the characteristic value (14. ip which decreases 
exponentially with the magnetic field, the inverse time (14. 7p grows and depends on B with 
the power law ~ B^. This increase comes from the e| factor in equation (14. 7p . which reflects 
the increased phase space available from the emission of phonons at high fields.— 



V. COMPARISON OF THE HYPERFINE COUPLING AND SPIN-ORBIT 
RELAXATION CHANNELS. DISCUSSION. 

Summing up the right-hand sides of Eqs. (I3.17p . (14. 2p . and (14. 6p . we find the total 
relaxation flux: 



dn^/dt = (n^ — n^^) 



^ 1 I srn I 



o I _ T-SrP c— ph 

^ V ^hf T^so Tso 



(5.1) 



As the inverse relaxation time is in fact proportional to n^, we characterize the relaxation 
process at a substantial initial excitation nx(0). The latter value experimentally is ~ 0.1 and 
under the assumed conditions T < 0.1 K and 5 > 10 T (where the equilibrium concentration 
ni^"* < 10"'') one finds the law nx(t) =nx(0)/[H-nx(0)t/2rtot], where 

Estimates of the r ' -values are possible if we specify material parameters included in 
formulae (l3.15p .( B?T]) and (14. 7p . Some of them have been already given in Sees. II-B and II- 
C and in Appendixes A and B. In addition we consider c = 5-10^ cm/s and =0.02555 meV. 
Other parameters related to modern wide quantum- well structures could be chosen as + 

= 10^^/5, , A = 50nm, A = 0.3 meV, and (i = 8.1nm (here B is assumed to be measured 



in Teslas; cf. also estimates in Ref. Il8l ). However, estimate of the effective spin-exciton 
mass Mx strongly depends on the finite thickness form-factor. There are experimental data 
where is found at comparatively low magnetic fields: (i) l/Mx~1.2meV at 5 = 2.27 T 
and = 1 in the 33 nm quantum well;— (ii) 1/Mx ~ 1.51 meV at i? = 2.69 T and u = 1 in 
the 23 nm quantum well;^^ and (iii) l/Mx~0.44meV at 5 = 2. 9T and z/ = l/3 in the 25 nm 
quantum well.— For these fields characterized by the inequality Ib > d, the 5-dependence 
should be approximately 1/Mx ~ B^^"^, but in the Ib < d strong field regime the inverse 
mass grows much more weakly with B. Based on these data, the semi-empirical analysis 
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using characteristic GaAs/AlGaAs form factors allows us to consider values 1/Mx~ 2meV 
at z/ = 1 and 1/Mx ~ 0.7 meV at z/ = 1/3 as the characteristic ones for the 10 T< B <25T 
range. (Note that at a given field B the estimate M~^|,,<i~z/'-M~^|,^i holds according to 
the semi-phenomenological theory.—) 

Numerical values of the characteristic inverse relaxation times are plotted in Fig.l. as 
a function of magnetic field. We remark that actual times should be longer by factor 
~ 2/r7,x(0) ~ 20 — 50 because of the non-exponential solution of equation fl5.2p . The in- 
dependence of the relaxation rate is non-monotonic. In the region 10 T < i? < 30 T the 
relaxation regime switches twice between the spin-orbit and hyperfine coupling dominance, 
taking maxima ^ 18T and ~ 12T in the u = 1 and z/ = 1/3 cases respectively. The rea- 
son that the hyperfine interaction becomes dominant is that for increasing magnetic field 
the nuclei remain disordered, while the random potential is effectively smoothed by the cy- 
clotron motion. At very high fields the spin-orbit interaction again dominates because of 
the increasing phase-space for the emission of phonons. On the basis of these estimates we 
conclude that the hyperfine coupling relaxation channel should be dominant approximately 
from 16 to 29 T in the u = 1 quantum Hall ferromagnet and from 11 to 24 T for u = 1/3. 
The latter case would seem to be more accessible to the experimental study of the hyperfine 
coupling relaxation mechanism, because usual electron concentrations in GaAs structures 
do not allow one to attain fields stronger 10 T in the u = l quantum Hall system. We note a 
feature of the hyperfine coupling relaxation: its rate is vanishing in the case of spin-polarized 
nuclei. This should distinguish the hyperfine coupling mechanism from that of spin-orbit and 
provide a test of the theory. If the nuclear spins could be fully polarized, then only spin-orbit 
relaxation would be important and there should be crossover between the regime limited by 
the random potential and the very high field regime of phonon emission. We emphasize also 
that our results should be valid in immediate vicinity of 1 or 1/3 fillings. Recent experi- 
ments show that if u differs by more than about 0.1 from these special values, one observes 
a two-mode spectrum of spin excitations - above and below the Zeeman gap.— Interaction 
of these two types of spin waves could considerably accelerate the relaxation. 

In conclusion, we have reported on a new spin relaxation mechanism in a spin polarized 
strongly correlated two-dimensional electron gas that appears at low temperatures and in 
strong magnetic fields. This mechanism is related only to the hyperfine coupling with GaAs 
nuclei and no other interactions are needed for this relaxation channel. The full calculation 
of relaxation displays a competition of the hyperfine coupling and spin-orbit relaxation pro- 
cesses, which can be summarized by equations ( I3.15p . (l4.ip . and (14. 7p . Under the assumed 
conditions the relaxation process occurs non-exponentially with time. The rate does not 
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depend on temperature but depends on the magnetic field non-monotonically as can be seen 
in Figure 1, which is plotted using estimated material and device parameters taken from 
experiment. The estimate of the hyperfine relaxation depends on the assumed randomness 
of the nuclear spins and a test of the theory would be to polarize the nuclear spins. 

S.D. thanks the Russian Fund of Basic Research and the LIA Condensed Matter and 
Theoretical Physics Program (ex ENS-Landau) for support, and the Laue Langevin Institute 
(Grenoble) for hospitality. 



Appendix A: Calculation of the hyperfine coupling parameters ylca and 

We proceed from formula An = (167r/iB/in/3/„)|M(i2„)p, — where /i„ is the nuclear 
magnetic moment, and u{Rn) is the conduction electron Bloch function at the nucleus. 
uiR) . a.u,ned to be normalized a= 1. where the mtegration . performed 

within the GaAs two atom unit cell having volume vq = 45.2 A'^. It seems to be the only 
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and subsequently cited by 



estimations of \u{Rqs)\'^ and \u{Ras)\'^ were done in Ref. 

other authors (cf. Ref. Isol ). Using these and the /i„ values for As and for the Ga®^ and 
Ga^^ stable isotopes: /ias = 1.44, fJ^Qges = 2.017 and /iQa'^i = 2.56 (in units of the nuclear 
magneton /i7v = 3.15 ■ IQ-^meV/G),^ we find ^^^69 ~ 0.038 meV, A^^n ~ 0.049 meV and 
Aas — 0.046 meV. Ratio of the Ga^^ and Ga^^ amounts in the semiconductor is considered 
to be equal to 3:2, therefore the result is 



J2 Al = 0.6 (Ag,69)^+ 0.4 {A^g.if+ (^^3-)^^ 4 ■ 10-3 meV^ . (A.l) 

within unit cell 



Appendix B: acoustic phonon life times T i and Tf 

If we take x, y, z to be the directions of the principal crystal axes, then for longitudinal 
phonons we obtain—"^ 

where 

" 2nhpc^' " I y 5hpc^ ^^"^^ 

(q and in this Appendix are considered to have common dimension.) Transverse phonons 
in a cubic crystal do not induce a deformation field.™ Actually, we need only the inverse 
time I/tj averaged over all directions of the transverse phonon polarization. If the transverse 
phonon distribution satisfies the condition that their polarizations are equiprobable, then for 
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either of the two polarization the averaging yields— 

We have used in Eqs. (B.2) and (B.3) common notations: H ~ 17.5 eV and eu ~ —0.16 C/m^ 
are the relevant deformation potential and piezoelectric constant of the GaAs crystal, p ^ 
5.3g/cm^ is the GaAs density, k, ~ 12.85 is the dielectric constant. As a result, we find 
Tf) ~ 0.8 ps and Tp ~ 35 ps. 



Appendix C 



The four-Q-operator expectation value is calculated with the help of Eq. (12. 4 p 
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where 



= (qi'x qi + q2'x q2)^ /2, and u' is considered to be equal i/ if z/ < 1, or 1 if 
the filling factor is integer. In the important case of integer u: {0\AqAl^\0) = 5q^o- Then, 
if calculating the matrix element (13. 2p . the (0|^...^^ |0) terms do not contribute to the 
probability transition (II. 4p owing to the energy conservation condition (13. ip . This means that 
the kinematic scattering would be determined only by the double-commutation expectation 
value (fin 



-^hf, Ql^ ) Qqj similar to the case of the dynamic scattering, cf. Eq. (13. 6p . 
In Eq. (13. 4p the action of the double-commutation term leads to the state 
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FIG. 1. Calculated inverse relaxation times as a function of magnetic field B from formulae ()3.15p . 
()4.ip and (j4.7p corresponding to hyperfine, l/rhf (solid line), spin-orbit with random potential 
I/t&^o^ (dash), and spin-orbit with phonon emission l/r|jf (dash-dot), respectively. Specific material 
parameters are given in the text. The bold solid line is the calculated combined inverse time (j5.2p . 



